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(54) PROCESS FOR CONSTRUCTING AND CODING LDPC CODE 

(57)Abstract: 

PROBLEM TO BE SOLVED: To easily construct an LDPC code for 
protecting a binary information string. 

SOLUTION: Respective information strings are composed of the N pieces 
of symbols decomposed into the N-M pieces of useful information symbols 
XI and the M pieces of redundant information symbols Ym and respective 
codes are defined by an inspection matrix A composed of N columns and 
M=N-K rows provided with the t pieces of non-zero symbols inside the 
respective columns. In this process, the same number of the non-zero 
symbols are allocated to all the rows of the inspection matrix A, the number i* — 
t of the symbols is an odd number as small as possible, the column is 
defined by a method that the optional two columns of the inspection matrix 
A are provided with only one non-zero value at maximum and the row is 
defined by the method that the two rows of the inspection matrix A are 
provided with only one non-zero common value. H 
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1.17 6 (2 H) 
1.16 7 (17 a— H) 

1. 143 (3 3-F) 

vhu^xA<^ff*5ctiK/^fctt^J 

[0 0 5 7] 
[#&] 



#^1 

for(x=0; x<n' - m + 1; x++) { 
pos[x] [0] =« O.- 
forCposCx] [l]«posIx] [01+1; pos[x] [1]<M-1; 
pos[x] [11++) { 

for(poa[xl (2]=pos[x] [1]+1; pes [x] [2l<M; 
pos[x] [2]++) { 

(if the conditions are not satisfied, 
continue 

otherwise, go to olc) 

} 

} 

(stop: impossible to find a suitable choice 
for pos [x] [0 — t-ll ) 



ok: ; 



} 
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(language C: the operator corresponds to EXCLUSIVE 

OR. 

/* Initialization of the b pass, of length M*/ 
for{i-M-k2; i<M; i++) 
b[i]=0; 

/* Calculation of N successive values of the 
impulse response of 1/A(z)*/ 
b[0] = 1; 

for(i=l; i<k2; i++) 

b[i] = bI(i+M-(k2-kl)) % M]''b[(i+M-k2) % M] ; 
for(i=k2; i<M; i++) 

b[i] = bli-{k2-kl)l^b[i-k2l; 

/* Arrange for there to be just one 1 

in the last k2 positions of b filtered by A(z)*/ 

weight = 0; /* all except 1 */ 

while (weight 1= 1) { 

/* Shift by one notch */ 

for(i=l; i<M; i++) 

b[i~l] = b[i]; 

b[M-l] = b[M-l-(k2-kl)]''b[M-l-k2l ; 

/* Verify */ 
weight = 0; 

for(i=M-k2; i<M; i++) { 

char sum = b [i] ''b [ (i+kl) %Ml ^b [ (i+k2) %M] ; 

if (sum) { 

shift « M - i; 
weight++ ; 

} 

} 

/* Particular case where M is forbidden */ 
if (weight ==0) 

return (M FORBIDDEN) ; 
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/* rightward final circular shift: 
b[il=b[(i - shift) % M] */ 
for{dec=0; dec<shift; dec++) { 

char tetnp = b[M-l] ; 

for(i=M-l; i>0; i--) 

blil=b[i - 1] ; 

b[0] = temp; 



return (OK) ; 
(language C) : 

/* Phase 1: calculation of M intermediate 
parities z. 

These parities are calculated by reading the 
successive columns of the coding matrix, namely 
At*l [Ml , . - .,A[*] [N] 

They are placed at the head of x ten^orarily */ 

ttdefine z x 

for(i=0; i<M; i++) 
z[il - 0; 

/* Loop over the n-m right submat rices of A*/ 

CO = M; 

cl a cO + P; 
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for(k = 1; k<=n - m; k++) { 
offset = 0; 
for(c=cO; c<cl? C++) { 
if (x[c] !=0) 

for(i=0; i<t; i++) { 

/* p ought to be 

offset + pos[i] - 

We decrement it by k2 to avoid 

shifting 

the array z before phase 3*/ 
p = offset + pos[k*t+i] - k2; 

if (p<0) 

z[p + M] = ztp + Ml^l; 

else 

if (p<M) 

z[p] = z[p]''l; 

else 

z[p - M]=z[p - Ml^l; 

} 

offset = offset + m; 

} 
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CO = cl; 

cl » cl+ P; 

} 

/* Phase 2: calculation of the last k2 parity 

symbols */ 

ixbO = M - 1 - lc2; 

/*!: initialization of the last lc2 elements of y 
temp[0...1c2-ll=y[M-l, M-2, ...M-k2l*/ 
for(k=0; k<k2; lc++) 
temp[kl = 0; 

1*2.: copy over the last k2 elements of b 
reglO. . .k2-l] = b [M-k2 . . .M-l] */ 
for(i=0/ i<k2; i++) 
reg[i] = finB[i] ; 

/*3: iterative calculation of the last k2 symbols 
*/ 

for(i=0; i<M; i++) { 

/* b(i] ={100 ...}"b[i-{k2-kl)pb[i-k2] 
with b[i-k2l . , .b[i-ll - reg [0 . . - t2-l] 

We must verify that: 

b[-k2] + b[kl-k2] + b[0] - 0 30 

b[-2] + b[kl-2] + b[k2-2] - 0 
b[-l] + b[kl-l] + b[k2-l3 = 0 
b[0] + b[kl] + b[k2l « 1 
b[l] + b[l+kll + b[l+k2] = 0 

... */ 

if (i"k2) 

input = 1; 
else 

input = 0 ; 
bi = input^reg [0] *reg [kl] ; 
for(k-l; k<k2; k++> 

reg[k - U = reg[k] ; 
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reg[k2 - 1] « bi; 

if(bi I- 0) 

for(lc=0; k<k2; k++) 

if(z[(ixbO - k + M) % M]!«0) 
tenq?[k] temp [k] ''I; 

ixbO = ixbO + 1; 
if (ixbO==*M) 
ixbO = 0; 

} 

/*4: The z values have already been left-shifted 
to avoid overwriting. Otherwise, it would be 
necessary to do: 



for(k=0; k<M - k2 ; k++) 
z[k] = zlk -I- k2] ; 



Copy over temp to the end of y */ 
#define y x 

for(k=0; k<k2; k++) 

y [M - 1 - kl = temp [k] ; 

/* Phase 3: calculation of y[M-k2-l, M-k2-2 , . . . , 0] 
y[k+k2-k2] + y[k+k2-kll + y(k+k2-0l + z[k+k2l = 0 
y [k] goes to x [k] 
z[k+k2] is in x[k] 
Hence : 

xlk+k2-k2] + x[k+k2-kl] + x[k+k2-0l + x[k] = 0 
i.e. : 

x(k+k2-k2] - - (x[k-i-k2-kll + x[k+k2l + x[kl) 
*/ 

for(k ^ M-k2-l; k>=0; k--) 

y[k] = y[k+k2-kl]^y[k+k2]''z[k] ; 

[03] 01coia^J<7)m2fe-9->^VhU^X<7)MXP 0T*^o 
^^-^ h U ^X'\<^^;l/-y^b^^T0Tfe'5o [0 8] 0 6 v=^y^^ V U ^7.^1^^ 

[04] c(D^m\zWi^'ya±7.(Dmi(D^mm\zm 0t*^o 

h U A^^f 0T^^o 9 ] 0 6 <hl^«|COX-r ^V^>^^ h U ^X^tin-T 

[0 5] z:(DmM\z^^':fa±x<Dm2<D^^m\zVt 50 m^&^o 
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